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ABSTRACT 

Aims. We analyze the first X-ray observations with XMM-Newton of 1RXS J070407.9+262501 and 1RXS 180340.0+401214, in 
order to characterize their broad-band temporal and spectral properties, also in the UV/optical domain, and to confirm them as 
Intermediate Polars. 

Methods. For both objects, we performed a timing analysis of the X-ray and UV/optical light curves to detect the white dwarf spin 
pulsations and study their energy dependence. For 1RXS 180340.0+401214 we also analyzed optical spectroscopic data to determine 
the orbital period. X-ray spectra were analyzed in the 0.2-10.0 keV range to characterize the emission properties of both sources. 
Results. We find that the X-ray light curves of both systems are energy dependent and are dominated, below 3-5 keV, by strong 
pulsations at the white dwarf rotational periods (480 s for 1RXS J070407.9+262501 and 1520.5 s for 1RXS 180340.0+401214). In 
1RXS 180340.0+401214 we also detect an X-ray beat variability at 1697 s which, together with our new optical spectroscopy, favours 
an orbital period of 4.4 hr that is longer than previously estimated. Both systems show complex spectra with a hard (up to 40 keV) 
optically thin and a soft (85-100 eV) optically thick components heavily absorbed by material partially covering the X-ray sources. 
Conclusions. Our observations confirm the two systems as Intermediate Polars and also add them as new members of the grow- 
ing group of 'soft' systems which show the presence of a soft X-ray blackbody component. Differences in the temperatures of the 
blackbodies are qualitatively explained in terms of reprocessing over different sizes of the white dwarf spot. We suggest that systems 
showing cooler soft X-ray blackbody components also possess white dwarfs irradiated by cyclotron radiation. 

Key words. Stars: binaries: close - Stars: individual: 1RXS J070407.9+262501, 1RXS 180340.0+401214 - Stars: novae, cataclysmic 
variables - X-rays: stars - Accretion, accretion disks 



1. Introduction 

Intermediate polars (IPs) are a subclass of magnetic cataclysmic 
variables (CVs) in which the rotational period P w of the white 
dwarf (WD) is not synchronized with the orbital period Pq of 
the binary system (P w <k Pq). The detection of a modulation at 
the WD spin period in the X-ray regime represents a signature 
of magnetic accretion. The material lost from the secondary star 
generally passes through a disc and is channelled onto the mag- 
netic polar regions of the WD where a strong shock develops 
(lAizulll973h . Hard X-ray radiation is emitted by the hot post- 
shock gas that cools via thermal bremsstrahlung while it set- 



* Based on observations obtained with XMM-Newton, an ESA sci- 
ence mission with instruments and contributions directly funded by 
ESA Member States and NASA, and with the Observatoire de Haute 
Provence (CNRS), France. 



ties onto the WD surface. Because of the highly asynchronous 
WD rotation and of the lack of detectable polarized emission 
in the optical and near-IR in most systems, IPs were believed 
to pos sess weakly magnetized WDs (B ~ 0.1 - 10 MG. IWarnert 
dl995h ). A distinct soft (kT ~ 20-50 eV) X-ray blackbody emis- 
sion was discovered by the ROSAT satellite in a handful of IPs 
dMason et alJll992t lHaberl et alJll994t lHaberl & Motor] 1 19951: 
Bur witz et al.ll 1996b . some of them also showing polarized op- 
tical radiation indicating field strengths up to 20-30 MG (e. g. 
iPiirola et al.l d 19931 120081) ). This spectral component resembles 
that observed in the X-ray spectra of Polars, the other subclass of 
magnetic CVs, that instead possess synchronously rotating and 
strongly magnetized (B ~ 10 - 230 MG) WDs. 

Our current understanding of the origin of soft X-ray emis- 
sion and energ y balance in Polars has gr eatly improved since the 
early model of Lamb & Masters (1979). It has become evident 
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that the soft and hard X-r ay emissions in at least t hese systems 
are largely decoupled (see lBeuerm ann 1999> l2004h . most of the 
soft X-rays arising from accretion of dense blobs ('blobby accre- 
tion') which penetrate deep in the WD atmosphere. Irradiation 
of the WD atmospheric polar regions by hard X-rays and cy- 
clotron radiati on emitted i n a sta nd-off shock has been recently 
discussed by Koni g et ail d2006l) . yelding to the conclusion that 
most of the reprocessed light appears in the far-UV rather than 
in the soft X-rays. 

The similarity with the Polars raised the question on 
whether the 'so ft IPs' are the 'true ' progenitors of low-field 
Polars (see also iNorton et al.l d2004l) ). However, very recently 
a soft and rather absorbe d X-ray emission is being found in 
a growing number of IPs ( Haberl et"aT1 20021; Ide Martino et alj 
2004, 2006a b, 2008; Evans & Hellierl2007HStaude et all2 008). 



Whether all IPs possess such a component and whether it bal- 
ances the primary emission from the post-shock flow is a new 
aspect that deserves investigation. 

Furthermore, in recent years the number o f candidate mem- 
bers of the IP group has gr own quite rapidly ( Wou dt & Warned 
120041: IGansicke et al.ll2005l) almost doubling the number of pre- 
viously known systems as well as enlarging the range of asyn- 
chronism. While this could s uggest an evolution towards syn- 
chronism (INorton etal.1 12004). the confirmation of new candi- 
dates is essential to understand the evolutionary link between 
IPs and Polars. 

1RXS J070407.9+262501 (hereafter RXJ0704) was listed 
in the ROSAT A ll-Sky Survey Source Catalogue (RASS, 
Voge s et al] {1999)) as a faint X-ray source, with a count rate of 
0.19cts s" 1 in the 0. 1-2.4 keV PSPC range and ROSAT hardness 
ratiosQHRl = -0.40 and HR2 = +0.55 that suggested the pres- 
ence of a soft component i n the X-ray spectrum. RXJ0704 was 
then identified as a CV by Gansi cke et al.l {2005) and proposed 
to be an IP from the detection of optical photometric and spec- 
troscopic periodicities at 480.71 s and ~ 250 min, respectively. 
The short period variability was ascribed to the spin period of the 
WD, while the other was identified as the orbital period of the 
binary system. The optical light curve folded at the putative WD 
rotational period was double peaked, suggesting the presence of 
two equally contributing accreting poles. However, no other X- 
ray observation was performed on this source in order to study 
its X-ray properties and hence to confirm its membership in the 
IP class. 

1RXS 180340.0+401214 (hereafter RXJ1803) was also dis- 
covered in the RASS as a faint (PSPC count rate 0.18 cts s~') 
source with a very soft X-ray spectrum (HR1 = -0.18 and 
H R2 = -0.02). Later, it was also identified as an IP candidate 
bv lGansicke et al.l (|2005) from its optical photometric variability 
at a 1520.5 s period that was interpreted as the rotational period 
of the accreting WD. A spectroscopic period of 160.2 min was 
inferred from the analysis of radial velocity variations in the Hq, 
emission line wings and identified as the binary orbital period, 
thus placing this CV in the 2-3 hr period gap. The two periods 
would then indicate a system with a weak degree of asynchro- 
nism. Also for this source no further X-ray observations were 
available. 

In this work we report XMM-Newton observations of 
RXJ0704 and RX1803 aiming at determining for the first 
time their X-ray variability as well as their X-ray broad-band 



1 ROSAT hardness ratios 1 and 2 are defined as HR1 = 
(ft40 ^:^'° 40) ^d HR2 = ^^fZ' °° ) . where ( A - B) is the 
raw c ount rate in the A- B energy range expressed in keV iMot ch et alj 
1996). 



spectral properties. For RXJ1803 we further complement the 
XMM-Newton data with optical spectroscopy acquired at the 
Observatoire de Haute Provence (OHP, France). 



2. The observations and data reduction 

2.1. The XMM-Newton observations of RXJ0704 

RXJ0704 was obs erved with the XMM-Newton observatory 
(Hansen et al.l l 2001l) on October 4 200 6 (OBSID: 0401650101) 
with the EPIC-pn (IStriideret al.ll200ll) and MOS dTurner et al.l 
1200 ll) cameras operating in full frame mode with the medium 
filter for total exposure tim es of 10.2 ks and 12.2 ks, respec- 
tively. The RGS instruments (Iden Herder et al .2001) were oper- 
ated in spectroscopy mode for an expo sure time of 13.9 k s. No 
useful data were collected with the OM (Ma son et al.l2 001). Due 
to a failure of the XMM-Newton mission operation centre sys- 
tem, the exposure was interrupted. RXJ0704 was then observed 
again on March 23, 2007 (OBSID: 0401650301) with the same 
instrument configurations. The OM was operated in fast imaging 
mode using the UVM2 filter that covers the ultraviolet spectral 
range 2000-2800 A. The total exposures were 8.1 ks for EPIC- 
pn, 9.6 ks for the MOS, 9.9 ks for the RGS and 6.2 ks for the 
OM. A summary of the observations is reported in Table Q] All 
the raw data have been reprocessed using the standard pipeline 
process included in the SAS 7.0 package. EPIC light curves and 
spectra were extracted using a circular region of 30"centred on 
the source. The same aperture radius was used to extract back- 
ground light curves and spectra on the same CCD chip of the 
target. To improve S/N, we selected up to double pixel events 
with zero quality flag for the EPIC-pn data, while for the EPIC- 
MOS cameras we increased the selection level up to quadruple 
pixel events. 

The EPIC data of October 2006 were characterized by a 
very low and rather constant background level (0.014 cts s and 
0.002 cts s -1 for EPIC-pn and MOS respectively). In order to 
obtain sufficient S/N ratios and use x 1 statistics in the spectral 
analysis, the average EPIC spectra were rebinned to obtain 20 
counts in each spectral bin. The photon redistribution matrix and 
the ancillary region file were created by using the SAS tasks rm- 
fgen and arfgen. The RGS1 and RGS2 first order spectra are too 
faint to allow a reliable analysis and, therefore, are not presented 
in this work. 

The March 2007 X-ray data were instead severely plagued 
by high background radiation for more than half of the duration 
of the observation. We then windowed the EPIC data to exclude 
those epochs of high (> 0.5 cts s~') background count rate, leav- 
ing us with the initial 2.9 ks of useful data, thus losing ~ 70% 
of the whole observation. The extracted EPIC light curves and 
average spectra were obtained using the same procedure as de- 
scribed before. Also in this case the RGS data did not provide an 
useful spectrum. The OM-UVM2 background subtracted light 
curve binned in 20 s intervals was obtained with the task om- 
fchain. The average count rate was 0.7 cts s _1 , corresponding to 
an instrumental magnitude UVM2 = 16.2 and to an average flux 

of 1.5 x 10~ 15 erg cirrV 1 A"' in the 2000-2800 A band. For a 
comparison, the continuum flux of the optical spectrum obtained 

bv IGansicke et al.l (120051) is ~ 2.2 x 10~ 15 erg cm" 2 s" 1 A~' at 
3800 A. 
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Object 


Instrument 


Date 


UT (start) 


Net exposure time (s) 


Net count rate (cts s ') 


RXJ0704 


EPIC-pn 


2006-10-04 


05:52 


10194 


1.036 ±0.010 




EPIC-MOS 




05:22 


12240 


0.232 ± 0.005 




RGS 




04:59 


13908 


0.032 ± 0.002 




EPIC-pn 


2007-03-23 


16:03 


2900 


1.463 ±0.023 




EPIC-MOS 




15:41 


4450 


0.364 ±0.011 




RGS 




15:40 


4721 


0.047 ± 0.001 




OM-UVM2 




15:49 


1239 


0.749 ± 0.025 








16:45 


1239 


0.822 ± 0.026 








17:11 


1240 


0.669 ± 0.023 








17:37 


1241 


0.634 ± 0.023 








18:03 


1241 


0.630 ± 0.023 


RXJ1803 


EPIC-pn 


2007-08-31 


02:18 


18457 


0.503 ± 0.005 




EPIC-MOS 




01:56 


19356 


0.148 ±0.003 




RGS 




01:56 


18737 


0.020 ± 0.001 




OM-B 




02:05 


1839 


1.645 ±0.030 








02:41 


1840 


1.808 ±0.031 








03:17 


1842 


1.850 ±0.032 








03:53 


1841 


2.040 ± 0.033 








04:29 


1839 


2.076 ± 0.034 




OM-UVM2 




05:05 


1840 


0.263 ± 0.012 








05:41 


1840 


0.239 ±0.011 








06:17 


1839 


0.222 ±0.011 








07:29 


1839 


0.255 ±0.012 



2.2. The XMM-Newton observations of RXJ1803 

RXJ1803 was observed on August 31, 2007 
(OBSID:0501230101) with all EPIC cameras operated in 
full frame mode and with the thin filter. The total exposure times 
were 20.0 ks and 21.7 ks for the EPIC-pn and MOS cameras, 
respectively. Of the two RGS instruments, only RGS1 could 
acquire data for 21.9 ks, while RGS2 was blocked in a 'set-up' 
mode during the observation. The OM was operated in fast 
imaging mode using sequentially the UVM2 and B filters, the 
latter covering the spectral range 3900^4900 A, for 9.2 ks each 
(see Table[T|i. The EPIC data were reprocessed and, then, filtered 
to avoid a flare event that occurred during the initial ~ 2000 s of 
the observation. Light curves and spectra were extracted using 
a radius of 30"centered on the source. We also applied similar 
pixel selection levels as done before for the EPIC-pn and MOS 
cameras and rebinned the spectra to obtain 25 counts per bin. 
Also for this source, the RGS 1 first order spectrum was too faint 
to allow useful analysis. Due to a failure in the omfchain routine 
of SAS 7.0, OM data were reprocessed at the MSSL using 
new, yet unreleased, routines. Unfortunately, there was still a 
processing problem in one of the exposures with the UVM2 
filter that cannot be solved. However, we obtained OM-B and 
OM-UVM2 background subtracted light curve binned in 10 s 
intervals. The average count rates were 1.88 cts s in the B 
band and 0.24 cts s in the UVM2 band. The corresponding 
instrumental magnitude were B = 18.6 and UVM2 = 17.3, 

while the average fluxes were 2.4 x 10~ 16 erg cm 4 s A and 
5.2 X 10~ 16 erg cirT 2 s _1 A , respectively. 

Heliocentric corrections were applied to the EPIC and OM 
data of both sources. 

2.3. OHP spectroscopy of RXJ1 803 

RXJ1803 was observed at the OHP during two nights on June 
29 and July 1, 2006 for 4.3 hr and 4.22 hr, respectively. 
Long slit spectra were obtained with the Carelec spectrograph 



dLemaitre et alj 1 19901) placed at the Cassegrain focus of the 
193 cm telescope and using a CCD EEV (2048 x 1024 pix- 
els) detector with 13.5 /mi pixel size. A 133 A/mm grating was 
used, with a slit width of 2", leading to a wavelength coverage 
of 3600-7200 A at a FWHM resolution of ~ 5.7 A. All exposure 
times were set at 1520 s, corresponding to the known optical pul- 
sation, to smear any effect from the spin variability. The obser- 
vations were performed in nearly photometric conditions with 
seeing typically ranging from 2" 5 to 3". We obtained 9 spec- 
tra of RXJ1803 each night. Standard reduction was performed 
with the ESO-MIDAS package, including cosmic rays removal, 
bias subtraction, flat-field correction and wavelength calibration. 
The wavelength calibration was checked on sky lines and all ra- 
dial velocity measurements have been corrected from the small 
instrumental shifts measured on the O I 5577 A line and from 
the Earth motion. Times have been converted in the heliocentric 
system. Flux calibration has been performed using the standard 
stars BD+28 421 1 and BD +33 2642 . 

3. Analysis and results 

3.1. RXJ0704 

We have analysed the X-ray and UV XMM-Newton data of 
RXJ0704 to search for variability at the optical period as well 
as to study its broad band spectral properties. 

3.1 .1 . The X-ray variability 

The light curves from all EPIC cameras were extracted in the 
range 0.2-15.0 keV with a 20 s time resolution. A clear signa- 
ture of a periodicity on a time scale of a few dozen of minutes 
is detected in both the October 2006 and March 2007 EPIC-pn 
data. Since EPIC-MOS light curves are very similar to those ob- 
tained with the EPIC-pn camera in both epochs, though with 
lower statistics, we only show the EPIC-pn data analysis. Fourier 
analysis was restricted only to the October 2006 data (Fig. [T]) be- 
cause of the short duration of the March 2007 observations. We 
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Fig. 1. Power spectra of the October 2006 EPIC-pn data of 
RXJ0704 in selected energy ranges. From bottom to top: 0.2- 
15.0 keV, 0.2-0.5 keV, 0.5-2.0 keV, 2.0-10.0 keV. The spin 
(oS) frequency, its harmonics and the orbital frequency (£2) are 
marked with vertical dotted lines. 



Fig. 2. Light curves of RXJ0704 folded at the 480 s spin period 
using the times of the dip quoted in the text. Top and central pan- 
els: full-band (0.2-15.0 keV) EPIC-pn light curves for October 
2006 and March 2007. Bottom panel: the OM-UVM2 spin pulsa- 
tion of March 2007 detrended from the low frequency variability 
(see text for details). 



found significant peaks at to ~ 180 day and at 2a>, but the 
strongest peaks occur at 3a> and Aid. There are also weak signals 
at 5a> and 8w, indicating a very structured X-ray light curve. A 
fit was performed on the October 2006 EPIC-pn and MOS light 
curves by using a composite function with up to six sinusoids to 
account for all the detected harmonics. We let the fundamental 
frequency a> free to vary, while the other five frequencies were 
constrained to be integer multiple values of at, i. e. 2a>, 3oj, 4u>, 
5co and 8w. In this way, we found a> = 179.88 + 0.06 day -1 for 
EPIC-pn, 179.7 ± 0.3 day" 1 for MOS1 and 179.85 ± 0.15 day" 1 
for MOS2. Because of the higher accuracy of the EPIC-pn data, 
we adopt P w = 480.3 + 0.2 s. This value is c onsistent within 
3cr wi th the optical photometric determination of Ga nsicke et alj 
(2005) and hence can be safely ascribed to the WD rotational 
period. 

Fig. |2] shows the full band EPIC-pn light curves of both 
epochs folded at the spin period. The light curves appear very 
complex and are dominated by two features: a dip, that we arbi- 
trarily placed at phase zero, and a bump at phase ~ 0. 15. The dip 
represents a fiducial marker of the spin pulsation. We hence de- 
rive the following times of the dip: HJD = 2454012.74140(7) for 
October 2006 and HJD = 2454183.17362(5) for March 2007. 
The spin pulse is more structured in March 2007, especially in 
the phase range 0.15 - 0.5. The dip and bump appear to be less 
defined, while an additional peak can be recognized at <p ~ 0.4. 
This latter could also be present in the October 2006 data, but 
with much lower amplitude. 

The Fourier analysis was also performed on the October 
2006 EPIC-pn light curves extracted with the same 20 s bin- 
ning time in selected energy bands: 0.2-0.5 keV, 0.5-2.0 keV 
and 2.0-10.0 keV (see Fig.Q]). While no sign of variability can be 
detected above 2.0 keV, significant peaks are evident in the soft 
band periodograms, where the spin frequency and its first three 



harmonics are clearly detected. The October 2006 folded light 
curves are shown in Fig. [3] The dip feature is clearly observed 
below 2 keV but the bump appears only below 0.5 keV. The 
March 2007 energy dependent light curves confirm the general 
behaviour observed in October 2006, with no sign of variability 
above 2 keV and a dip whose depth increases at lower energies. 
The fractional depth (relative to the average count rate level) of 
the dip is, within errors, slightly larger in October 2006 than in 
March 2007 (75 + 7% and 66 + 12%, respectively). Similarly 
the bump centered at phase ~ 0.15 has a fractional amplitude of 
87 + 14% in October 2006, while its amplitude is 62 + 10% in 
March 2007. Both the dip and the bump are extremely localized 
in phase (A0 ~ 0.1), therefore the regions responsible for these 
features should have a very small extension. 

In the lower panel of Fig. [3] we also plot the hardness ratio 
defined as the ratio between the 0.5-2.0 keV and 0.2-0.5 keV 
count rates (the hardness ratio obtained with the other energy 
bands does not provide useful information). The count rate ratio 
of the soft bands is double-humped showing a hardening during 
the dip as well as at <p ~ 0.3 - 0.5 that corresponds to the 'flat' 
part of the spin pulse. This clearly reflects the presence of the 
soft bump and the dip. 

3.1.2. The UV light curve 

The March 2007 OM-UVM2 light curve shows hints of a quasi- 
sinusoidal modulation on the order of few hours. However, this 
cannot be quantified because of the short duration of the obser- 
vation. To better study the short term variability we removed the 
long term variations from the UV light curve using a third or- 
der polynomial and folded the residual at the X-ray spin period 
with the epoch of the dip found for March 2007. As shown in the 
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Fig. 3. RXJ0704 data of October 2006. Top panels: EPIC-pn 
light curves in selected energy ranges folded at the 480 s spin 
period using the time of the dip quoted in the text. Bottom panel: 
EPIC-pn hardness ratio in the soft bands. 



lower panel of Fig. [2] the pulsation is quasi-sinusoidal. A sinu- 
soidal fit gives the maximum at phase 0.22+0.03, hence partially 
overlapping the X-ray bump and the additional peak at <p ~ 0.4 
observed at that epoch, and a modulation fraction of 12 + 5% 
(defined as 2A// ave , where A is the amplitude of the fitted sinu- 
soid and / ave is the average count rate). We also recall that the 
white light optical pulsation observed bv lGansicke et al. (2005) 
was instead found to be double-humped. 

3.1.3. Spectral properties of RXJ0704 

The October 2006 (Fig. @> and March 2007 average EPIC-pn 
and combined MOS spectra were analyzed with the XSPEC 12 
package in the range 0.3-10.0 keV. A simple model consisting 
of an absorbed optically thin plasma MEKAL with solar abun- 
dances does not adequately fit the spectra. For October 2006 
data, a zero-width Gaussian centered at 6.4 keV is also required, 
although the fit is still poor (xl = 8.8). Both October 2006 and 
March 2007 spectra show an excess of counts below 1 keV and 
the addition of a blackbody (IcTbb ~ 80 eV) drastically improves 
the fits (xl = 0.97 and 0.70, respectively). For the October 
2006 spectrum, the temperature of the optically thin plasma is 
unconstrained (&7mek > 44 keV), while for the March 2007 
spectrum it is A:7mek = 1 IkeV. However, a multi-temperature 
plasma, expected to be present in the post-shock region, does 
not improve the fit quality. If the abundance is left free to vary, 
it assumes an unrealistic value Az = 1 .7 for the October 2006 
data, while for the March 2007 spectrum we obtain Az = 0.6, 
but still within errors compatible with solar abundances. We 
therefore fix Az to the solar value. The fits further improve by 
including an additional dense absorber covering the source by 
~ 40%, as detailed in Table [2] The inclusion of a denser ab- 
sorber is justified by a F-test at significance level of 99.9%. 
The low column density of the total absorber, when compared 



Fig. 4. The October 2006 EPIC-pn spectrum of RXJ0704 is 
shown with the best fit model discussed in the text and parame- 
ters reported in Table[2](for the sake of simplicity, we do not plot 
the combined MOS spectrum). Each spectral component that 
contributes to the best-fit model is shown separately, for clar- 
ity. The bottom panel shows the residuals expressed in terms of 



with the total galactic absorp tion in the direction of the source 
(A^gai = 8.15 x 10 20 cirr 2 . [Dickey & Lockmanl dl990h ). sug- 
gests a likely interstellar origin. On the other hand, the dense 
absorber should be confined within the binary system and likely 
due to pre-shock material as suggested by the energy depen- 
dence of the spin light curve. No major differences are found 
in the spectral parameters between the October 2006 and March 
2007 data (see Table except for the lack of detectable iron 
6.4 keV line in the March 2007 low S/N spectrum and a signif- 
icant change in the normalization of the blackbody component, 
which increases by a factor of two in March 2007. This change 
could explain the increase in the source flux observed between 

< 10~ 12 erg crrr 2 s and 
t 



the two epochs (F°* mi koV 

pMar. 

r 0.2-10.0 keV 



= 5.3 x 10 12 erg cm 2 s respectively). 



We also performed a phase-resolved analysis of the October 
2006 EPIC-pn spectra (the inclusion of the MOS spectra does 
not greatly improve the statistics) extracted around the very lo- 
calized dip (</> = 0.9 - 1 .05) and bump (0 = 0.1 - 0.2), as well as 
in the flat part (</> = 0.3-0.8). Since the X-rays are essentially un- 
pulsed above 2 keV, we used the best-fit model of October 2006 
keeping fixed the MEKAL parameters to the values found for 
the average spectrum. We also fixed the column density of to- 
tal absorber as this is not expected to vary along the spin cy- 
cle. The Gaussian component used to account for the iron line 
at 6.4 keV is not required to fit the spectra extracted arount the 
dip and the bump because of the poor S/N ratio. Unfortunately, 
with our data many spectral parameters are badly constrained, as 
shown in Table [3] Anyhow, we find a substantial change in the 
covering fraction of the partial absorber between the flat portion 
of the pulse and both the dip and the bump. The values of the 
normalization of the blackbody component found in the flat part 
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Table 2. Spectral parameters of the best- fit model to the EPIC-pn and combined MOS average spectra of RXJ0704. Errors indicate 
the 90% confidence level of the relevant parameter. Metal abundance is kept fixed at the solar value (see text). 



Epoch 


Total Absorber 
Nul 

(10 20 cm" 2 ) 


Partial Absorber 

Nut Cov - Frac.^ 
(10~ cm" 2 ) 


Blackbody 
kT BB N BB t 
(eV) (10- 5 ) 


MEKAL 

£7mek A'mek'I 
(keV) (10~ 3 ) 


Iron Line 

EW f 

(eV)" 


^(^/d.o.f.) 


Oct. 2006 
Mar. 2007 


1 o +u - 6 

i.u _ 05 

0.5 ±0.4 


1 9 +ly 
1 2 +L2 


33+ ulj 

u,JJ -0.0! 

54+ - 07 

-0.09 


84 ±3 
88 ±5 


3 +oy 
6 3 +L7 


> 44 
11+' 3 


2.3 ±0.3 
2 6 +0 - 4 


141 ±73 


0.87 (360/415) 
0.63 (121/193) 



Column density of the total absorber. 
Column density of the partial absorber. 
Covering fraction of the partial absorber. 
Normalization constant of the BBODY model. 
Normalization constant of the MEKAL model. 
Equivalent width of the 6.4 keV iron line. 



Table 3. Spectral parameters obtained from fitting the October 
2006 phase-resolved EPIC-pn spectra of RXJ0704 extracted in 
the phase intervals quoted in the text with the model shown in 
Table 



1RXS J1803+4012 



Parameters 


Dip 


Bump 


Flat 




F 0.2- 10.0 keV 


2.71 


5.15 


4.13 




(10~ 12 erg cm" 2 s" 1 ) 










Nu (10 23 cm" 2 ) 


2 l +u 


1 3 +LI< 

-0.7 


2 7 +1J 

^■'-0.7 




Cov. Frac. 


0.73 ± 0.04 


0.48^ 


0.44 ± 0.20 


m 


kT BB (eV) 


84 +u 


82 ±4 


85 ±2 




N BB (10- 5 ) 


T ^> + l.o 

-'"'-0.8 


7 7 +L3 

'•'-1.0 


3.6 ± 0.2 






xi Cr z /d.o.f.) 


1.15 (26/23) 


1.03 (63/61) 


1.00(297/298) 





and in the dip are consistent, within the errors, but there is an 
evident increase in the bump. Hence, it is conceivable that the 
spin pulse of RXJ0704 is due to an increase in the local absorp- 
tion (dip feature) and to an increase in the projected area of the 
blackbody component (bump). 



3.2. RXJ1803 

Similarly to RXJ0704, the XMM-Newton data were also an- 
alyzed to search for periodicities in the X-ray, UV and opti- 
cal ranges. The spectroscopy acquired at OHP was analyzed to 
search for variability in the radial velocities of Balmer lines. 

3.2.1. The X-ray periodicities in RXJ1803 

The EPIC-pn light curve (EPIC-MOS light curves were too 
noisy to provide useful information) was extracted in the full 
energy range of the camera with a resolution of 20 s and Fourier 
analyzed. The power spectra in the 0.2-15.0 keV range, reported 
in Fig. [5] show the presence of a strong peak at 56.6 day - and 
an additional weaker peak at 49.8 day -1 . A sinusoidal fit com- 
posed with one single frequency gives then a period of 1528+6 s. 
This is fairly consistent with the hig hly accurate opt i cal p eriod 
of 1520.510 + 0.066 s determined bv lGansicke et all (120051) . al- 
lowing us to confirm it as the true spin period of the WD. The 
residual light curve when fitted with a second sinusoid gives a 
period of 1697 + 22 s. If this frequency represents the sideband 
lo - CI (commonly named as the beat), using either the X-ray 
period or the optical period we derive an orbital frequency of 
5.64 + 0.69 day 1 or 5.91 + 0.66 day -1 , respectively. This gives 
an orbital period of 4.25 + 0.52 hr or 4.06 + 0.45 hr which does 
not match either the orbital period of 2.67 hr determined by 
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Fig. 5. Power spectra of the EPIC-pn light curve of RXJ1803 in 
selected energy ranges. From bottom to top: 0.2-15.0 keV, 0.2- 
0.5 keV, 0.5-1.0 keV, 1.0-3.0 keV and 3.0-10.0 keV. The spin 
(of), and the orbital sideband (at - Q) as determined by optical 
photometry and our spectroscopy (see text) are marked with ver- 
tical dotted lines. 



iGansicke et al.l (120051) . or 2Pq. As it will be shown below, the 
orbital period is likely 4.4 hr, hence longer than previously de- 
termined. The detection of X-ray variability at the spin and beat 
frequen cies (with different amplitude s) is also observed in other 
IPs (see lNorton etaLlll992alfbf Il997l) . In RXJ1803, the ratio of 
power between spin and beat frequencies is ~ 3.8, implying that 
the beat modulation has an amplitude about twice lower than that 
of the spin. 

Light curves of the EPIC-pn were also extracted in different 
energy ranges with a binning time of 20 s and Fourier analyzed, 
as also shown in Fig. [5] The spin variability is clearly detected 
below 3 keV, as well as the beat frequency. The ratio of power 
between spin and beat peaks is ~ 4.8 in the 0.2-0.5 keV, ~ 3.6 in 
the 0.5-1. OkeV and ~ 2.8 in the 1. 0-3.0 keV, suggesting that the 
variabilities at the two frequencies do not have the same energy 
dependence. Because of the higher accuracy of the optical period 
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Fig. 6. The spin pulse of RXJ1803 in selected energy ranges of 
the EPIC-pn and in the B and UVM2 bands. All the light curves 
are folded at the 1520.5 s optical spin period using the time of 
maximum quoted in the text. 



1RXS J1803+4012 




0.2 - 



1 2 3 4 5 

Time (hr) 

Fig. 7. Analysis of the X-ray spin pulse of RXJ1803 in the range 
0.2-3.0 keV. Upper panel: the time dependence of the phase of 
the spin maximum. Lower panel: the time dependence of spin 
amplitudes. 



determination, we folded the spin light curves at the 1520.5 s 
period (Fig. [6} using the time of maximum obtained from a si- 
nusoidal fit with fixed frequency: HJD = 2454343.6124(4). The 
folded light curves are quasi-sinusoidal with a dip at <p ~ 0.8 
which is more evident in the softest 0.2-0.5 keV band. The pulse 



amplitude (see Sect. [3~L2l is 72 ± 8% in the 0.2-0.5 keV range, 
57 + 7% in the 0.5-1 .0 keV and 46 ± 6% in the 1 .0-3.0 keV band, 
indicating a strong energy dependence especially below 1 keV. 
The hardness ratios in the three bands however do not substan- 
tially vary along the spin cycle, indicating that the spin pulse is 
dominated by aspect angle variations. 

The beat variability suggests that the shape of the spin pul- 
sation is variable along the binary period. The XMM-Newton 
observation spans about 5.6 hr, likely covering more than one 
orbital period. We then extracted the spin light curve in six in- 
tervals each sampling ~ 3000 s (i. e. about two spin cycles). 
While amplitudes are rather constant within errors, the spin pulse 
changes phase with time (see Fig. [7j». Altough we cannot estab- 
lish whether this variability is periodic, the presence of the beat 
variability suggests it might be the case. 

3.2.2. The UV and optical light of RXJ1803 

The OM-UVM2 light curve does not reveal any clear trend, 
while that in the B filter shows both a long term and short term 
variability. Count rates steadily increase by a factor of ~ 1.4 
during the 2.6 hr of observation with the B filter. This could re- 
flect an orbital dependence, suggesting, as found in the X-ray 
data, that the orbital p eriod is longer than previously estimated 
(Gan sicke et al.l 120051) . We then removed the long term trend 
in the B-band light curve using a third order polynomial. The 
folded spin light curves in the two bands are shown in Fig. [6] 
Although very noisy, they appear to be different: in the UV no 
modulation is apparent within errors, while the B band light 
curve is quasi-sinusoidal, with an amplitude of 3 + 10%. The B 
band p ulse is single peaked, as also detected by iGansicke et al.l 
d2005l) . and is broadly consistent in phase with the X-ray pulse. 

3.2.3. The X-ray spectrum of RXJ1803 

The EPIC-pn (Fig. [8) and combined MOS spectra averaged 
over the whole observation were fitted using a model consti- 
tuted by an absorbed MEKAL with solar abundances. This sim- 
ple model does not satisfactorily fit the spectrum (xl = 2.2). 
A slight improvement is achieved by using a multi-temperature 
plasma CEMEKL (xl = 1.95) with kT m!lx ~ 25 keV and power- 
law slope a fixed at unity. In both cases, there is an excess of 
counts below 0.5 keV and above 5 keV, but the hydrogen col- 
umn density of the total absorber drops to zero, thus indicating 
that this component is not required. As a comparison, the to- 
tal galactic absorption in the direction of RXJ1803 is Nu, S a\ = 
3.4x 10 20 cirT 2 dDickev & Lockm anl 19901) . Leaving metal abun- 
dance free, the fit still remains unsatisfactorily high (xl = 1 .4) 
giving sub-solar abundances Az = 0.2. A great improvement 
(xl = 0.91) of the fit is found with a composite model (Model A 
in Table |H consisting of a blackbody and the multi-temperature 
plasma, which decreases the excess of counts at low energies. 
The blackbody temperature is found at IcTbb ~ 100 eV. A fur- 
ther improvement is found by adding a partial covering (40%) 
dense (N H = 2.8 x 10 23 cirT 2 ) absorber (Model B in Table @). 
The presence of this component is significant at the 98.6% level 
with the F-test. Here we note that the metal abundance of the 
optically thin plasma is still sub-solar (Az = 0.4) and that the 
slope of the power-law emissivity function was fixed to unity, 
because when left free it assumed the value a = 0.67 but consis- 
tent with unity within errors. A similar quality fit (xl = 0.89, 
Model C) is obtained by using two MEKAL components in- 
stead of the CEMEKL, supporting the idea that the optically thin 
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Table 4. Spectral parameters of the models used to fit the EPIC-pn and combined MOS average spectra of RXJ1803, as described 
in the text. Errors indicate the 90% confidence level of the relevant parameter. 





Partial Absorber 

Nut Cov - Frac-l 
(10 S cm" 2 ) 


Blackbody 
(eV) (10- 6 ) 


A zL 


Optically thin pi 
kTi Nil 
(keV) (10" 3 ) 


asma^ 
kT~± 

(keV) 


N 2 h 
(10~) 


Xt Ct 2 /d.o.f.) 


A 
B 

C 


3 +3 0.4 ± 0.2 

5 +5 5 +03 

J_ 2 "- J -0.2 


98 ± 8 2.3 + 0.2 

QC + 10 Q Q + 1.4 

yj_ q J - J -o.6 
104 ± 10 4.2^ 5 


j 2+0.4 


"•^-0.2 

2 +0 - 2 


1.0 ±0.2 


3 +L0 


40+i" 

21! 

12™ 


1 61 +,,uy 
i.ui_„ 07 

T a +0-8 

1 6 +26 

i ' u -0.7 


0.91 (448/495) 
0.89 (437/493) 
0.89 (436/491) 



" The optically thin plasma is represented by a CEMEKL (Models A and B) with fixed slope a = 1 or by two MEKALs (Model C). 

b Column density of the partial absorber. 

c Covering fraction of the partial absorber. 

d Normalization constant of the BBODY model. 

e Metal abundance. 

' Normalization constant of the first MEKAL model (only for Model C). 

s In Models A and B this is the maximum temperature of CEMEKL model. 

* Normalization constant of the CEMEKL (Models A and B) or the second MEKAL model (Model C). 



plasma is not isothermal. The parameters of the partial absorber 
and the blackbody are consistent with that of Model B, but badly 
constrained. The temperatures of the two MEKAL components 
are 1 keV and 12 keV, the latter being consistent within errors 
with the maximum temperature of the multi-temperature plasma 
found in Model B. The low temperature plasma is instead re- 
quired to fit a shallow bump located aroud 1 keV and ascribed to 
a complex of spectral lines likely dominated by a strong O VIII 
line at 19 A. Unfortunately, we cannot inspect the RGS because 
of the lack of data from RGS2. The observed X-ray flux in the 
0.2-10.0 keV range is F Qt2 _io.o keV = 1-84 x 10~ 12 erg ctrT 2 s -1 . 
A comparison with the observation carried out during RASS 
shows that RXJ1803 was about 1.5 times brighter in the soft 
ROSAT band during the XMM-Newton observation. In 2007 the 
source is also less soft than when detected during RASS, since 
the ROSAT hardness ratios would have been HR1 = -0.09 and 
HR2 = +0.12, suggesting that the soft X-ray component was less 
prominent in the XMM-Newton observation than during RASS. 

We also performed a phase-resolved analysis of the EPIC-pn 
spectrum extracted around the pulse maximum (<p = 0.85 -1.15) 
and minimum (<p = 0.35-0.65), in order to investigate the mech- 
anism responsible for the pulsation. Model B shown in Table [4] 
was used, keeping the metal abundance fixed to the value ob- 
tained for the average spectrum and fixing the power-law index 
or = 1. The parameters of the CEMEKL component are very 
similar in both fits and are completely consistent with the value 
found in the average spectrum. This was expected since we do 
not observe any variabilty above 3 keV. Given quoted errors, 
the main variations arise in the normalization of the blackbody 
component, which increases at the spin maximum by a factor 
of 2.7 with only a weak evidence that this occurs in the par- 
tial absorber column density. The lack of a significant variabil- 
ity in the hardness ratios suggests that the X-ray spin pulse in 
RXJ1803 is likely dominated by changes in the projected area 
of the blackbody component. It is interesting to note that the ra- 
tio between the spin and the beat amplitudes is larger at lower 
energies, which could suggest a much lower contribution of the 
reprocessed soft X-ray component to the beat variability. 

3.2.4. The spectroscopic period of RXJ1803 

The mean optical spect rum of RXJ1803 was found similar to 
that already reported by Gan sicke et alj d2005l) and is typical of 
magnetic CVs with strong emission lines of the Balmer series, 
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Fig. 8. The EPIC-pn spectrum of RXJ1803 fitted with Model B 
discussed in the text (see Table [4j. The combined MOS spec- 
trum, also used in the fit, is not shown. Each spectral component 
that contributes to the best-fit model is shown separately. The 
bottom panel reports the residuals expressed in terms of cr. 

Table 5. Best-fit spectral parameters of the phase-resolved 
EPIC-pn spectra of RXJ1803 extracted in the phase intervals 
quoted in the text. Model B shown in Table|4]has been used. 



Parameters 


Maximum 


Minimum 


f 0.2-10.0 kcV 


1.94 


1.55 


(10~ 12 erg cm" 2 s" 1 ) 






N H (10 23 cm" 2 ) 


2 9 +4 3 


j 4 +4.2 


Cov. Frac. 


0.5 ± 0.2 


6 + ^ : ' 
106!" 

1 6 +u 
'•"-as 


kT BB (eV) 


103 + J8 
4 3 +3 - 


A'bb (10- 6 ) 


kTmax 




13 + f 


NcEM GO" 3 ) 


3 4 +1 ° 


3 2 +1 - 9 


Xi (X 1 / d.o.f.) 


1.27 (131/110) 


0.80 (50/62) 
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Fig. 9. Power spectrum of the H a radial velocities. The best value Fig. 10. Radial velocities of the H a line folded with the orbital 
is found at Pq ~ 5.45 day -1 (~ 4.40 hr). Secondary maxima are period Pq = 4.402 hr. The best sine fit is also shown (line). 
2 day aliases. 

The XMM beat value of 1697 + 22 s is therefore consistent with 
Table 6. Parameters of the orbital modulation of the emission the above orbital period, 
line radial velocities. Errors are at the lcr confidence level. 



Lines 


y (km s 1 ) 


K(km s~') 


Phase " 




-94.6 + 4.9 


69.2 + 6.8 


0.0 (fixed) 


H, 


-70.3 + 6.4 


72.2 + 8.7 


0.02 + 0.13 


Hy 


-39.3 + 9.1 


37.9+12.6 


-0.01 ±0.15 



" Phase of the blue-to-red zero crossing 



Hell (4686 A) and He I (4471, 5875, 6678 and 7065 A), su- 
perimposed on a relatively blue continuum. The radial veloci- 
ties for the main lines were measured using a single Gaussian 
least square fit procedure and periodicities were searched using 
a Fourier analysis. The power spectrum for the strongest H ff line 
gives the most accurate period determination (see Fig. 0. The 
maximum power is seen at 5.45 + 0.27 day -1 with secondary 
peaks which are 1 and 2 days aliases, corresponding to the night 
separations. The characteristics of the radial velocity line modu- 
lation were further determined from ax 2 sine fit. The minimum 
X 2 corresponds to a period of Pn = 4.402 + 0.014 hr, where 
the error bar is at a lcr level, computed for two independent 
parameters. The orbital ephemeris for RXJ1803 is determined 
as T n = HJD 2453917.361(20) + 0.18342(57) E, where T n is 
the predicted heliocentric time of the blue-to-red radial velocity 
transition. 

The radial velocities of the H„ line, folded at the above best 
orbital ephemeris, are shown in Fig. [10] and the parameters of 
the corresponding best sine fits for the strongest emission lines 
are given in Table [6] Though the 4.40 hr period is the most sig- 
nificant, we cannot exclude the 2 day aliases at 4.029 + 0.014 hr 
and 4.846 + 0.014 hr. These observations therefore do not con- 
fir m the period gap value 2.670 ± 0.008 hr previously reported 
by IGansicke et al.ld2005l) . We note however that a large scat- 
ter is visible in their Fig. 22 showing the radial velocities of 
Hq, wings folded at the 2.67 hr period and that these are only 
marginally represented by a sinusoidal fit at their period. A sig- 
nificant distortion of the line profiles may have been introduced 
by the 600 s and 1800 s exposures used, not commensurate with 
the spin value. Using the 4.40 hr value for the orbital period and 
the optical spin value (1520.51 s), the u-Q. and o> + O beat peri- 
ods are predicted at 1690.7+8.4 s and 1393.4+6.9 s, respectively. 



4. Discussion 

The X-ray analyses of RXJ0704 and RXJ1803 reveal 
strong pulses a t the optical periods previously identified by 
IGansicke et al.l d2005l) that now can be interpreted as the true 
rotational periods of the accreting WDs and hence confirming 
these 2 CVs as IPs. This increases the current roster of con- 
firmed IPs to 3 1 systems Q. Both sources show pulsations that 
are strongly energy dependent in the soft X-ray portion of the 
spectrum. RXJ0704 and RXJ1803 also show common spectral 
properties characterized by a hard and a soft emission compo- 
nent. The latter property adds these two systems to the small but 
growing group of 'soft IPs'. Now we discuss their pulse prop- 
erties and spectral characteristics highlighting their common as- 
pects. 

4.1. RXJ0704 

The 480 s spin pulsation observed below 2 keV is highly struc- 
tured with a strong dip feature and a bump, while no modu- 
lation is detectable at higher energies. The spin pulse is very 
similar to that o bserved in PQ Gem (de Mart ino et al.l 2004: 
lEvans et ai1l2006l) . which also shows a marked dip that preceeds 
the spin maximum. The dip in both sources is due to an in- 
crease of absorption from an intervening column along the line 
of sight. While for PQ Gem the geometry of the system an d 
its accretion flow is rather well determined dPotter et alJI 19971) . 
that in RXJ0704 has still to be established. Notwithstanding 
this, the strong similarity of the pulse may suggest that also in 
RXJ0704 the dip could result from the obscuration of the mag- 
netic pole produced by the accretion flow along field lines (cur- 
tain, iRssenetal] (11988)) which preceeds the pole itself. In this 
configuration, the material couples at large distances and trav- 
els out of the orbital plane producing large absorption effects. 
At later phases, when the curtain moves away from the line of 
sight, the X-rays are less absorbed and the WD area of the ac- 
creting pole is still at its maximum visibility. This is also cor- 
roborated by the variability seen at UV wavelengths, where the 

2 For a comprehensive list of IPs see also 
http : / /asd . gsf c . nasa . gov/Ko j i . Mukai/ iphome/iphome . html 
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maximum is observed in correspondence with the X-ray maxi- 
mum. RXJ0704 has been observe d in a previous epoc h to show 
a double-sine optical modulation dGansicke et al.ll2005l) . thus in- 
dicating that both poles are active at some epochs. To have both 
of them come into view the magnetic dipole should be relatively 
inclined (m < i, where m is the magnetic colatitude and i is the 
binary inclination). Eclipses are not observed in RXJ0704, im- 
plying i $ 70°. Furthermore, the lack of modulation in the hard 
X-rays suggests that either the accretion shock does not pass be- 
hind the limb of the WD or there are two poles equally contribut- 
ing to the hard X-rays. This last is favoured in the framework of 
two visible poles. Wether material couples at large radii, pro- 
ducing very localized absorbing structures, need to be explored 
with further observations such as phase resolved spectroscopy 
and multi-band photometry. 

The spectral analysis confirms the presence of a blackbody- 
like soft X-ray component which strongly dominates the pulse 
below 0.5 keV. This component is relatively hot and should 
originate in a relatively small area of the WD polar cap. In 
fact, the normalization constants of the blackbody in both the 
epochs of observation imply an emission area that is in the range 
4 - 7 x 10 12 dl m cm 2 E|. The bolometric flux of the blackbody 

is Feb = 2.6 - 5.6 x 10" 12 erg crrr 2 s -1 . On the other hand, 
the optically thin emission has a bolometric flux 4.6 -1.4 times 
larger than that of the soft component in the two observations. 
Therefore, to obtain an estimate of the mass accretion rate, we 
need to include both components. Here we note that the UV flux 
is modulated at the spin period by 12% and, when dereddened 
by Eb-v = 0.015 (as de rived from the hyd rogen column den- 
sity of the total absorber. iRvter et ail d 19751) ). is about three or- 
ders of magnitude higher than the extrapolated blackbody flux 
in that range. Hence, the accretion luminosity is estimated as 
L accr > L B b + ihard + iuv = 1.1 x 10 33 af 2 00pc erg s" 1 . The up- 
per temperature value of the optically thin component of 24 keV 
provides a lower limit to the shock temperature, which in turn 
would give a WD mass of O.52M0. We then assume a canonical 
WD mass value of 0.6M© WD and derive an accretion rate of 
~ 10 16 c/g 00pc g s~' . At a distance of 800 pc it is lower by one or- 
der of magnitude than the 2.6 x 10 17 g s _1 expected for its orbital 
pe riod JwTrnerll 19951). 

iNorton et alj d!999l) proposed that fast rotators like 
V709 Cas (Pu = 313 s), YY Dra (P u = 529 s) and 
V405 Aur (P m = 545 s) possess weakly magnetized WDs. 
RXJ0704, with its short 480 s spin period, joins these fast ro- 
tators. Its magnetic moment could be estimated assuming that 
the WD is accreting from a disc. Th e strong spin variabil- 
ity indicates that accretion is disc-fed dNortoniri993l) . Hence, 
the condition for disc formation and truncation at the magne- 
tospheric radius implies that R mdg ~ R co , where the corota- 
tion radius is the radius at which the magnetic field rotates 
with the same Keplerian frequency of the inner edge ot the ac- 
cretion disc. These are defined as R co = (GM^Pl ) /4n z \ ' 

and# mag = 5.5 x 10 8 (M WD /M©) 1/7 R~ 2n L^ n ^ cm, where 
Rg is the WD radius in units of 10 9 cm, L33 is the lumi- 
nosity in units of 10 33 ergs -1 , and fi^o is the WD magnetic 
moment in units of 10 30 G cm 3 . For a O.6M0 WD, we find 

3 For RX J0704 and RXJ 1 803 we adopt lower limits to their distances 
of 800 pc and 500 pc, respectively. These are obtained using the 2MASS 
K-band magnitudes, adopting for the secondaries a spectral type M3 ± 2 
for RXJ0704 and M3± 1 for RX J 1803, that are predicted for their orbital 
periods of 4. 17 hr and of 4.4 hr (Knigae 2006), and assuming that these 
donors totally contribute to the K-band flux. 



H ~ 1.0 x 10 32 o?800pc G cm 3 . The spin-to-orbit period ratio of 
RXJ0704 is Paj/Pn ~ 0.032, thus it is a relatively highly asyn- 
chronous system which poses the question of whether this sys- 
tem is sp inning at equilibriu m. Spin equilibria have been com- 
puted by Nort on et al.1 d2004) for systems with a fixed mass ratio 
q = 0.5 and accreting at the secular value for their orbital pe- 
riods. For Poi/Pa = 0.03 and an orbital period of 4.2 hr, it is 
predicted that equilibrium is reached for a magnetic moment of 
/ieq ~ 5 x 10 32 G cm 3 . While this could suggest that RXJ0704 
is not spinning at equilibrium, this issue remains open until the 
distance and the mass ratio remain undetermined. 

4.2. RXJ1803 

The X-ray timing analysis of RXJ 1803 has revealed the pres- 
ence of a strong pulsation, with an amplitude of 54%, at the pre- 
viously identified optical period, which confirms this system as 
an IP. We also detect a weaker modulation (by a factor of ~ 2) 
at a period of 1697 + 22 s that we interpret as the beat period. 
This is consistent with the value obtained using the optical pho- 
tometric spin period and our new spectroscopic orbital period, 
i. e. 1681.88 ± 0.56 s. Hence, both our new spectroscopy and 
X-ray data show tha t the binary period is much larger than pre- 
viously determined (Gansicke et al. 2005). With the new 4.4 hr 
orbital period, RXJ 1803 is outside the orbital period gap and its 
spin-to-orbit period ratio lowers to P^/Pn ~ 0. 096, which is 
well w ithin the range found for most systems (see INorton et al.l 
12004 . 

X-ray pulsations at the beat frequency occur when the ac- 
cretion stream is directly channelled to the magnetic poles 
of th e WD, without passing through the disc dWvnn & King! 
119921) . However, the detection of both spin and beat period- 
icites suggests an hybrid accretion mode, i. e. occurrin g si- 
multaneously via a disc and disc-overflow dHellier| [l995). The 
proportion between the two accretion mechanisms can be in- 
ferred from the relative amplitudes of the corresponding X-ray 
variabilities. Disc-overflo w was also theoretically predicted by 
lArmitage & Livid dl998l) . To date, only the IP V2400 Oph is 
known to be a pure stream-fed accretor. A few other systems 
(such as FO Aqr, AO Psc, TX Col) have shown both spin and 
beat variabilities, with different amplitudes. RXJ 1803 thus joins 
these few systems with disc accretion representing ~ 66% of the 
total flow. 

The X-ray spin pulse below 3 keV is single peaked and 
energy dependent, while no detectable modulation is found at 
higher energies. Both hardness ratios in the soft bands and spec- 
tral analysis show that the spin pulse is likely dominated by as- 
pect variations rather than photo-electric absorption. We also 
find changes of the phase of pulse maximum on a timescale 
of hours, commensurable with the orbital period. The X-ray 
beat variability does not appear to possess the same spectral 
dependence of the spin, although it is also observed below 
3 keV. The spectrum of RXJ 1803 is well described by an op- 
tically thin multi-temperature plasma with maximum temper- 
ature kT m . dx ~ 30 - 40 keV which can be regarded as the 
shock temperature. A blackbody component at 95 eV also con- 
tributes to the soft portion of the spectrum by 27% of the to- 
tal observed flux below 1 keV. Both components are absorbed 
by a dense (A^h = 3 x 10 23 cm -2 ) local material covering 
~ 40% of the X-ray emitting region. If the latter is mostly re- 
lated to the overflowing material, this could explain the differ- 
ent energy dependence of the beat and spin variabilities. The 
blackbody emission is also relatively hot, with a very small 
area on the WD pole (~ 1.0 x 10 11 t/ 2 00 c cm 2 ). Its bolomet- 
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ric flux is F BB = 2.9 x 10~ 13 ergcrrT 2 s _1 which results to 
be 14 times smaller than that of the optically thin component. 
To evaluate the accretion luminosity we include both hard and 
soft X-ray components, as well as the modulated flux in the 
UV/optical bands. However, the UV range does not show a 
modulation. We therefore include only the B band modulated 
flux of 7.3 x 10~ 14 ergcrrT 2 s _1 , corresponding to a luminos- 
ity 2.2 x 10 30 ^5 ()0pc ergs -1 . The accretion luminosity is then 

iacci- ^ ihard + L BB + L B = 1.3 x 10 32 af 2 00pc erg s -1 . To esti- 
mate the mass accretion rate we use a WD mass inferred from 
the shock temperature of 30-40 keV: M W d ~ 0.61 - 0.76 Mq. 
This gives an accretion rate M t 0.8 - 1 x 10 15 ^g 00pc g s _1 . If the 
source is indeed at 500 pc, this value is low and much smaller 
than the secula r value predict ed for its 4.4 hr orbital period, i. e. 
3 x 10 17 g s -1 dWarnerj|1995l) . Although this is a lower limit to 
the accretion rate and unless most of the accretion luminosity is 
not emitted at high energies, RXJ1803 seems to be a soft IP with 
a low accretion rate. 

The ratio between the spin and orbital period is 0.096, which 
is very close to the typical value of 0. 1 for IPs and suggests that 
RXJ1803 is spinning close to equilibrium. The dominance of 
the spin pulsation also indicates that RXJ1803 accretes predom- 
inantly via a disc. As done in Sect. 14.11 the condition for disc 
truncation at 7? mag gives a lower limit for the magnetic moment 
li ~ 1.3 x 10 32 fl?5oopc Gem 3 u sing the M and the WD mass 
derived for this object. Following Nor ton et al.l {2004), the mag- 
netic moment predicted for a system in spin equilibrium with 
Pu/Pn = 0.096 and P n = 4.4 hr is fx = 2 x 10 33 G cm 3 . Further 
observations to constrain the multicolour variability as well as 
the distance of this system are needed to draw firm conclusions 
on the mass accretion rate and the magnetic moment of this sys- 
tem. 

4.3. Soft X-ray IPs: an emerging new class? 

Our XMM-Newton observations of two new IPs allowed us 
to detect a soft X-ray blackbody component in these sys- 
tems. The ROSAT satellite was the first to discover three IPs 
with a significant soft X-ra y componen t, nam e ly PQ Gem, 
V405 Aur and UU C o l (Mason et all fl99l faaberl et al l 
Il994t iMotch et alj 119961; iBurwitz etalj 1 1996b . Such a com- 
ponent was also recognized in a few other systems ob - 
served with BeppoSAX (V2400 Oph, Ide Martino etaT] ( 12004 ) 
and with th e high sensitive XMM-Newton data (NY Lup, 



lHaberletal.1 d2002h : 1RXS J2133.7+5107, Ide Martino etal] 



(2006b); 1RXS J173021. 5-055933. Ide Martino et al.l (|2008) 



1RXS J062518.2+733433.IStaudeetal.M2008h). A systematic 
search for this component in the XMM-Newton archival data 
has further added three systems EX Hya, GK Per and WX Pyx 
dEvans & Hellierl 12007). The current roster of soft X-ray IPs 
amounts to thirteen systems, which is now a statistically sig- 
nificative sample, comprising 42% of all the confirmed IPs. 
Therefore, 'soft IPs' seem to represent a well defined class 
among IPs. This makes the soft X-ray emission not a sole 
characteristics of Polar systems which deserves an explanation. 
Moreover, the properties of the soft X-ray component in IPs still 
have to be defined. With thirteen systems, this is now possi- 
ble. We compiled Table [7] a listing of the blackbody temper- 
atures and softness ratios (defined as the bolometric flux ra- 
tio F so ft/4Fhard, see lRamsav & Cropper! d2004 ) as derived from 
spectral analyses available in literature. Despite the great spread 
due to different models applied and/or different spectral data, 
there is clear evidence that the soft X-ray blackbody in IPs cov- 



ers a wide range of temperatures (30-100 eV), that is much 
wider than in the Polar systems (see also Fig.fTTTi. In most cases 
the temperatures are high and larger than those found in the 
Polars with only a few systems falling in the same range. These 
hot blackbodies could arise from small cores of the accretion 
spots o nto the WD poles. As already noted in Ide Martino et all 
(2008), although the emitting areas are subject of great uncer- 
tainties, mostly due to the distance, systems showing hot black- 
bodies like NY Lup, 1RXS J173021.5-055933, RXJ0704 and 
RXJ1803 have fractional areas / ;S 10~ 5 - 10~ 6 which are two- 
three orders of magnitude lower than those found in IP s with 
cooler blackbodies such as PQ Gem (Jam es et al. 1 12002b . since 
the inferred blackbody temperature is an average value over the 
heated spot, the wide range of temperatures might reflect accre- 
tion spots of different sizes. 

Strong evidence is found in most Polars that the soft X-ray 
component arises in the WD photosp here primarly due to the 
heating by accretion of dense blo bs (Kuijpers & Pringiel [l982; 
IWoelk & Beuermanrifl996T: lBeuermannlll999ll2004 . with most 
of the reprocessed primary (hard X-rays and cyclotro n) radia- 
tions mainly emitted in the far-UV dKonig et al.ll2006b . As seen 
in Fig. QT| the softness ratio of IPs is generally smaller than that 
observ ed in the Polars and, as already noted by IE vans & Hellierl 
(l2007h . this indicates that blobby accretion does not occur in 
IPs. We further stress that the presence of a disc in these sys- 
tems indeed makes this accretion mode very unlikely. Most IPs 
also have softness ratios below 0.25, lower than that predicted by 
the standard model L, epr « 0.5 (Lb re m + ^cyc) dLamb & Mastersl 
1979), suggesting that the soft X-ray emission does not fully ac- 
count for the whole reprocessing. 

lEvans & Hellierl d2007b propose that the appearance of this 
component is primarly due to geometric factors because the 
cool material in the accretion curtains heavily absorbes the soft 
X-rays when viewing along the curtain. They note that those 
IPs with light curves dominated by absorption dips tend not to 
show a blackbody component. While the high density local ab- 
sorbing material certainly adds further complications in IPs this 
might not be the whole story. Our finding of a wide range of 
temperatures, and the likely different sizes of the heated spots, 
could be explained if bremsstrahlung and cyclotron radiation 
both play a role in irradiating the WD surface. T he fact that 



(Shakhovskoi & Kolesnikovl 11997 


), V2400 Oph dBucklev et al.l 


1 1995b and 1RXS J2133.7+5107 ( 


Kataiainen et al.l 2007b, were 



found to show optical/near-IR circular polarization, suggests that 
cyclotron radiation cannot be neglected. Here we note that it can- 
not be excluded that other 'soft IPs' are also circularly polarized 
but have not been searched in polarized light yet. 

Beamed cyclotron radiation originates in the higher parts 
of the post-shock region and, thus, can heat large spots. 
Bremsstrahlung is instead irradiated isotropically from the dense 
post-shock regions closer to the WD surface and hence irradi- 
ates smaller areas. Therefore we might expect that for higher 
magnetic field IPs, cyclotron cooling is not totally negligible 
and therefore these systems would have larger accretion spots 
with respect to other IPs with lower magnetic fields. This in turn 
would imply that higher field IPs are expected to show on aver- 
age lower blackbody temperatures than the low field IPs. A test 
to this hypothesis would then require a systematic search for po- 
larized radiation in the soft IPs and, in particular, over a wide 
range of wavelengths from visual to IR. 

The d etailed and recent study on reprocessed radiation in 
Polars by Kon ig et all J2006) relies on the clear evidence of a 
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Table 7. Blackbody temperatures and softness ratios (see 
text) of 13 'soft IPs', as derived from spectral fits avail- 
able in literature. For some objects different values are re- 
ported as different models were applied and/or different 
spectra l data were used. Refer e nces a re (1) Ide Martino et al 




(4) lHaberl et al 



Evan s & Hellie 
T 



3 (120071). (3) Ide Martino et al . 
2h. (5) Ide Martino et alj d2006b). 



(6) Ide Martino et all 



H002l 

20081) . m lStaude et alj d2008l) . 
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0.009 
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V405 Aur 
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65 
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0.048 
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82 


0.02 
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MU Cam (2005) 
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MU Cam (2006) 


54 
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7 


RXJ0704 (Oct. 2006) 
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0.053 


This work 


RXJ0704 (Mar. 2007) 


88 


0.172 


This work 


RXJ1803 


95 


0.019 


This work 



large irradiated spot on the WD surface in AM Her that mainly 
emits in the far-UV, implying that most of the soft X-rays do 
not enter in the energy balance. For IPs, there is still a great 
uncertainty in whether the UV modulated flux comes from the 
irradiated WD pole or from the inner regions o f the curtain or 
both (e . g.lde Martino et al l dl999ll2006 a.'2008): Eisenbart et al] 
d2002l) : lEvans et all (120061) ). Detailed far-UV/UV coverage to 
identify the heated WD surface in IPs would be highly desirable 
to assess the energy balance in these systems. 



5. Conclusions 

XMM-Newton observation of RXJ0704 and RXJ1803 confirm 
these CVs as true members of the IP class. Both systems are 
dominated by strong pulsations at the WD rotational periods 
which dominate the soft portion of their spectrum, i.e. below 3- 
5 keV. With P w = 480 s, RXJ0704 joins the group of fast rotators 
and possess a relatively high degree of asynchronism {P^/Pn — 
0.04). RXJ1803 is, instead, a slow rotator P w = 1520.5 s. From 
our new optical spectroscopy and the X-ray data of this system 
we infer an orbital period of 4.4 hr, that is longer than previously 
determined. This gives a spin-to-orbit period ratio of ~ 0.1, typ- 
ical for IPs. 

The detection of strong spin pulses indicate that accre- 
tion occurs preferentially via a disc in both systems. However, 
while RXJ0704 only shows X-ray variability at the spin period, 
RXJ1803 also displays a weak variability at the orbital sideband 
{to - Q), implying that accretion also occurs via disc -overflow. 

Their X-ray spectra are compatible with composite emission 
consisting of an optically thin component (~ 1 1 -25 keV), which 
in RXJ1803 is better described by a multi-temperature plasma, 
reaching values up to 30^40 keV, as well as by blackbody emis- 
sion with a relatively high temperature (85-100 eV) heavily ab- 
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Fig. 11. The softness ratio of Polars (filled triangles) and of 'soft 
IPs' (other symbols) v ersus blackbody temperatu re. Values for 
Polars are taken from Ram say & Crop per (2004), while those 
for IPs from Table □ 



sorbed by dense material partially covering the X-ray source. 
This adds these two systems to the growing group of 'soft IPs', 
which indicates that the presence of a soft blackbody compo- 
nent is not solely a characteristic of the Polars. However, a ma- 
jor difference is in the temperatures of the blackbodies in the IPs 
which cover a wider range than that observed in the Polar sys- 
tems. We qualitatively try to explain the differences in terms of 
reprocessing over different sizes of the WD spot. We suggest that 
IPs showing cooler soft X-ray blackbody components also pos- 
sess WDs irradiated by cyclotron radiation. A further systematic 
observational approach addressing the detection of polarized op- 
tical/IR emission in IPs is highly desiderable. 
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